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Abstract 

It  is  shown  that  the  variational  method  of  dealing  with  the 
integral  equations  of  scattering  problems  is  equivalent  to  solving  the 
integral  equation  directly  by  Galerkin's  method  and  using  the  standard 
fonnula  for  the  amplitude  of  the  scattered  wave.  The  second  method  also 
satisfies  the  reciprocity  theorem.  It  is  therefore  suggested  that  the 
reciprocity  theorem  be  \ised  as  the  basis  of  approximation  without  the 
introduction  of  variational  formulas. 

The  error  involved  in  using  an  approximate  solvition  is  discussed 
and  it  is  shown  that  only  a  special  set  of  approximations  can  lead  to  ac- 
curacy at  low  frequencies.  Some  ways  in  which  bounds  for  the  error  may  be 
obtained  in  special  problems  are  also  given. 


Table  of  Contents 

1.  Introduction  1 

2.  The  equivalence  of  the  variational  method  2 
and  Galerkin's  method 

3.  The  accuracy  of  the  approximation  7 
U.  The  low  frequency  behavior  10 
5»  Boxinds  for  the  error  12 

References  18 


NEW  YORK  UNIVERSITY  _^ 

INSTITUTE  OF  MATHEMATICAL  SCIENCES 
LIBRARY 
25  WaveHy  Place,  New  Jfprk  S,  H,  Yj 


-  1  - 


1,  Introduction 

The  fact  that  an  integral  equation  can  be  formulated  as  a  variational 
principle  was  stated  as  early  as  188U  by  VolterraL  -■  but  it  has  not  been  used 
for  practical  calculation  in  scattering  problems  until  recent  years  when  it 
was  discovered  by  Schwinger  that  the  amplitude  of  the  scattered  field  is  close- 
ly related  to  the  quantity  whose  variation  has  to  be  considered.  (In  connec- 

r2] 

tion  with  small-amplitude  sound  waves,  see  Levine  and  Schwinger  ••  '',  and  in  con- 
nection  with  atomic  collisions  see  Schwinger »s  results  as  quoted  by  Kohn"^-*).  The 
essence  of  the  variational  method  is  that,  if  a  good  approximation  for  the  field 
is  inserted  in  the  variational  expression,  an  improved  approximation  for  the  scat- 
tered amplitude  should  result.  One  particular  way  of  choosing  the  approximation, 
that  given  by  Levine  and  Schwinger  '■  -* ,  is  to  expand  the  field  in  a  set  of  functions 
and  solve  the  simultaneous  linear  equations  resulting  from  the  variational  principle 
for  the  coefficients.  It  will  be  shown  in  the  following  that  this  is  exactly  equi- 
valent to  Galerkin's  method  of  solving  the  integral  equation  and,  furthermore, 
that  only  by  expanding  the  field  in  a  special  set  of  functions  can  accuracy  be 
obtained.  Also  tiiis  solution  satisfies  the  reciprocity  theorem. 

Another  way  of  approximation  is  to  choose  for  the  field  an  expression 
which  is  physically  plausible  and/or  mathematically  simple  and  insert  it  immed- 
iately in  the  variational  expression.  It  will  be  shown  that  this  is  equivalent 
to  choosing  an  approximation  T*iich  satisfies  the  reciprocity  theorem. 

It  appears,  therefore,  that  fundamentally  the  use  of  the  variational  method 
is  to  ensure  that  the  reciprocity  theorem  is  complied  with.  Since  the  reciprocity 
theorem  is  easily  written  down,  it  is  obviously  simpler  to  formulate  our  analysis 
in  terms  of  it  than  to  introduce  the  variational  principle.  The  advantage  is 
mainly  one  of  clarity  for,  once  the  equation (s)  to  determine  the  approximation 
have  been  obtained,  the  labor  required  to  calculate  the  scattered  amplitude  is 
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the  sane  by  either  method, 

2»  The  equivalence  of  the  variational  method  and  Galerkin's  method' 

The  mathematical  problem  encountered  in  scattering  may  be  stated  as 
follows:  given  a  function  f  and  operator  L,  find  a  function  g  such  that 

(1)  f  -  l€  . 

The  function  f  is  determined  by  the  incident  field  and  g  can  be  regarded 
as  a  suitable  source  distribution.  The  equation  is  supposed  to  have  a  unique 
solution.  Examples  of  equation  (l)  are 

f(r')  -  I   g(r)  T(r,  r')  dr 

for  one  of  the  boundary  conditions  in  diffraction  by  an  aperture,  and 


f (r  )  -  g(r  )  + 

J 


r 

V(r)g(r)Y(r.r  )dr 


S 
for  atomic  collisions  or  scattering  by  a  dielectric  obstacle.  In  these  equations 


T(r,  r')  -  ^ 


-iklr-r' 


Un|r-r  | 

S  is  *".ne  area  of  the  aperture,  T  is  the  volume  of  the  scatterer  and  7  is  the 
potential  of  the  scattering  particle,   (For  the  dielectric  obstacle  V  is  essen- 
tially  a  measure  of  the  dielectric  constant.) 

Often  it  is  the  scattered  field  at  some  distance  from  the  obstacle  and 
not  g  itself  which  is  of  primary  concern.  The  amplitude  A  of  the  scattered 
field  may  be  written 


I  Some  results  similar  to  those  of  this  section  have  been  obtained  independently 
by  Ilppmann'-  -J . 


« 


Eq,  (1)  and  the  following  analysis  are  applicable  to  vector  field  problems 
on  regarding  f  and  g  as  vectors  and  L  as  a  vector  operator  but  no  specific 
example  will  be  quoted. 
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(2)  A  -  (f',  g)  -  (g,  f') 

where  f  is  some  suitable  incident  field  and  (f  ,  g)  means  the  inner  product 
of  f  and  g.  For  example,  the  amplitude  of  the  far  field  in  the  direction  of 
the  unit  vector  n  is  obtained  correctly  if  we  define  (f  ,  g)  by 

(f',g)  -  j   f'(r)g(r)dr 
for  the  aperture  problem  or  by 

(f',  g)  -  I  f'(r)V(r)g(r)dr 


for  the  atomic  collision  and  take  f  ■  e  -  -  .  (A  different  choice  of  f 
would  lead  to  some  other  quantity,  the  near  field  for  example,  but  would  not 
invalidate  the  following  analysis  provided  that  it  was  a  possible  incident 
field. ) 

Let  g  be  such  that 

(3)  f  -  l€  . 

Assume  that  L  is  a  symmetric  operator  such  that 

(U)  (h^,  liig)     •     (hg,  l^i)  . 

Then,  from  (3),    (U)  and  (1), 

^^^  (f',  g)     -     (!«',  g) 

-  (g',   l€) 

-  (g',  f) 

A  -   (f  ,   g)      -      (f,   g    )     -      (g  ,    l€)    . 

This  is  the  reciprocity  theorem  which  exists  for  most  scattering  problems. 
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By  an  obvious  use  of  this  theorem,  we  can  write 

(6)  A  -  (^'>6)(^i6')   . 

(g,  !€*) 

It  can  be  confirmed  that  the  necessary  and  sufficient  condition  for  (1)  and 
(3)  to  be  satisfied  is  that  this  expression  for  A  is  stationary  for  independent 
small  variations  of  g  and  g  about  their  correct  values.  Other  variational 
principles  are  available  (see,  for  example,  Kohn^  -■  )  but  will  not  be  considered 
here. 

To  demonstrate  that  the  choice  of  an  approximation  to  make  (6)  station- 
ary is  equivalent  to  Galerkin's  solution  of  (1)  we  shall  consider  first  the 
special  case  in  which  f  ■  f  so  that  g  ■  g  •  This  corresponds,  for  example, 
to  the  problem  of  the  circular  aperture  in  a  normally  incident  plane  sound  wave 
with  the  scattering  observed  in  the  direction  of  propsigation  of  the  incident 
wave.  In  this  special  case 

(7)  A  .  ^^^     . 

(g,l€) 

Let  0  ,  G.,...  be  a  set  of  functions  each  of  which  is  independent  of 
the  others  and  each  of  which  satisfies  any  conditions,  e.g.,  edge  conditions,  im- 
posed on  g.  If  an  infinite  number  of  G  ,  G  ...  is  used,  the  set  is  also  required 
to  be  complete.  Now  try  an  approximation  G  to  g  given  by 

(8)  G  .  n  c^G^ 

n 

where  the  sum  may  involve  either  a  finite  or  infinite  number  of  terms  and  the 
constants  c  are  to  determined.  We  shall  now  prove  the  following 
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Theorem  1  The  necessary  and  sufficient  conditions  for  the  express- 
ion in  (7)  to  be  stationary  for  small  variations  of  G  are 

(9)         (^»  V  '  (°n'  ^^  (n-0,1,...). 

G  then  satisfies  (f,G)  =  (G,  LG)  and  A  -  (f,  G). 

The  importance  of  this  theorem  lies  in  the  fact  that  equations  (9) 

are  precisely  those  obtained  by  substituting  the  expansion  (8)  in  (1)  and 

taking  the  inner  product  of  both  sides  with  G  ,  which  is  Galerkin's  process. 

Then  A  may  be  calculated  directly  from  its  original  definition  (2)  because 

the  approximation  satisfies  the  reciprocity  theorem.  Obviously  this  direct 

method  of  setting  up  the  equations  is  to  be  preferred  to  determining  them  via 

the  variation  principle.  It  should  be  remarked  that,  if  we  make  a  one -terra 

approximation  so  that  G  »  c  G  ,  the  reciprocity  theorem  is  satisfied  if 

00 

(f,  G^)  «  c^(G^,  LG^)  and  (f,  G)  becomes  (f,  Gq)^/(Gq,  UG  )  which  is  the  same 
as  would  be  calculated  from  (7). 

We  proceed  now  to  prove  the  theorem.  The  first  variation  of  (7)  with 
g  ■  G  when  c  becomes  c  +  5  c  is 

(10)        -L 1 __L  . 

(G,  mf 


Suppose  now  that  equations  (9)  are  satisfied.  Then  (10)  vanishes  if 

(f,  G)  -  (G,  LG).  Now  (f,g)  -  (i,  ZI  <^nO   "  (^>  21  ^a^a)    °"  account 

of  (9);  hence 

(11)  (f,  G)  -  (G,  10) 

so  that  the  first  variation  is  zero.  The  equations  (9)  are  thus  sufficient 
for  (7)  to  be  stationary. 
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Suppose  now  that  the  first  variation  is  zero  for  independent  6c. 
Then 

(G,  LG)(f,  G^)  -   (f,G)(LG,  G^)         (n  -  0,  1,..,). 

These  equations  are  homogeneous  in  the  constants  c  ,  c^  . « .  .  Hence  if  there 

is  one  solution  there  are  any  nuinber  of  the  form  Co  ,  Cc,  •••  vrtiere  C  is  a 

o'   1 

constant©  There  is  therefore  a  solution  in  which 

(12)  (G,  LG)  -  (f,  G)  . 

Such  a  solution  satisfies  equations  (9)  and  also  the  reciprocity  theorem.  The 
necessity  for  the  equations  has  therefore  been  shown. 

On  accoimt  of  (11)  and  (12),  A  »  (f,G)  and  the  theorem  has  been  proved. 

The  corresponding  theorem  for  the  more  general  expression  (6)  will  now 
be  stated.  No  proof  will  be  given  since  it  runs  parallel  to  that  of  Theorem  1. 
We  take  G  to  be  the  approximation  to  g  and  try 

G  =  T^  c  G  , 

^—   mm' 
m 

this  series  containing  as  many  terms  as  (8).  Then 

Theorem  2  The  necessary  and  STifficient  conditions  for  the  express- 

ion  (6)  with  g  »  G  and  g  ■  G  to  be  stationary  for  small  independent 

t 
variations  of  G  and  G  are 

(f,  G*)  -  (IC,  G^)  (ra  -  0,  1,...) 

and 

(f',  G„)  -   (IT)',  OJ   . 

G  and  G  satisfy  the  reciprocity  relation  (f  ,  G)  -  (G,  IC  )  ■  (f,  G  ). 


The  value  of  A  is  either  (f,  G  )  or  (f  ,  G), 
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3»  The  accuracy  of  the  approximation 

To  determine  a  good  approximation  it  is  necessary  to  specify  precisely 
what  is  meant  by  good*  This  involves  stipulating  the  behavior  of  the  operators 

and  functions  in  some  manner.  We  shall  choose  the  behavior  which  is  appropriate 

f 

to  a  perturbation.  Attention  will  be  concentrated  on  the  case  f  ■  f  so  that 

A  ■  (f>  g)  and  the  reciprocity  theoran  takes  the  form  (f,  g)  »  (Ig,  g). 

Let  a  be  a  parameter  such  that 

(i)  for  some  range  of  a 

2 
L     5     L+aL.+a     L^+,,, 

f     -     f^  +  a  f     +  a^  f     +  ... 

O  1  2 

g     -     g^  +  a  g-L  ♦  a     gg  ♦  ... 
G     -     0°  ♦  a  G-^  +  a^  G^  ♦  .,.       , 

(11)    the  only  solution  of  L  h  -  0  is  h  ■  0, 

(iii)   (h^,  Lp  h^)  -  (Lp  h^,  h^)  for  all  p. 

When  a  is  small  we  are  writing  the  operator  L  as  one  which  is  slightly 
perturbated  from  L  .  If  a  measures  the  ratio  of  the  linear  dimensions  of  the 
obstacle  and  the  wavelength , and  L  is  the  static  operator, then  we  are  consider- 
ing  the  low-frequency  behavior.  It  has  been  shown  by  Magnus^-'  that  assumption 
(i)  is  justified  for  the  problem  of  a  low -frequency  plane  wave  falling  on  a 
circular  aperture. 

If  we  substitute  the  formulas  of  (i)  in  (l)  and  equate  the  coefficients 
of  powers  of  a  the  equations  to  determine  g  are 

(^^  ll  Vr  ^  "  ^P  (p-0,  1,...). 

r"0 
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If  the  inverse  of  the  operator  L  is  known^g  can  be  found,  in  principle,  correct 
to  any  power  of  a.  Assumption  (ii)  ensures  a  imique  solution* 

The  approximation  G  is  to  be  determined  so  that  the  reciprocity  theorem 
is  satisfied,  i.e.,  so  that  (G,  Ifi)  ■  (f,G),  Equating  the  coefficients  of  the 
like  powers  of  a  in  the  expansions  of  both  sides,  we  obtain 


or 


(lU) 


f:  F    (g^  l      q"  -  l      g ") 


on  account  of  (13  )• 

Now  we  are  going  to  approximate  (f,g)  by  (f,G).  Both  quantities  can  be 
expanded  in  powers  of  a  and  the  coefficients  of  a^  are  the  same  if 

''''         £,  (vs' °°- ^=.)  ■  °- 

Consequently  the  first  P  terms  of  (f,g)  are  given  correctly  by  the  approximation 

if  (1$)  holds  f or  all  p  <  P  -  1. 

Let  us  now  suppose  that  G  -  g  (m  <  M),  Then,  in  (lU),  we  can  replace 

m    ~ 

G^  by  g  in  all  the  inner  products  in  which  p  -r  <  M  for  in  these  the  second 
member  is  identically  zero.  Also  G^  may  be  replaced  by  g^  in  all  inner  products 
in  >rtiich  r  <  M,  By  this  means  G^  is  replaced  by  g^  in  all  inner  products  in  (lU) 
provided  that  p  <  2M  +1.  Hence,  after  a  use  of  assumption  (iii). 
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t     £'    (Vr-s  V  "'  -is)     ■    °  (P  S  2M  »  1) 

r-0  s"0   V*^  y 

or 

on  interchanging  the  order  of  summation*  This  equation  becomes,  on  taking 
account  of  (,13), 


£.   i^n-s'  ^'-^sJ   ■  °  (P<2M*1) 


which  is  the  same  as  (l$).  Thus  if  G°  •  g^^^  (m  <  M)  the  first  2M  ♦  2  terms 
of  (f,g)  are  given  correctly  by  the  approximation* 

If  we  carry  out  a  similar  analysis  with  p  ■  2M  +  2  we  find 


so  that  the  (2M+3)rd  term  of  (f,g)  is  given  correctly  only  if 


(16) 


(f''  -  %.i.  V"*'  -  V«.i)   ■  0  • 


This  result  enables  us  to  obtain  some  information  which  is  particularly  applicable 
to  the  low-frequency  behavior.  At  low  frequencies,  where  L  is  the  static  opera- 
tor and  a  is  purely  imaginary,  the  following  additional  assumptions  hold  in  most 
scattering  problems i 

(iv)  L  operating  on  a  real  function  produces  a  real  function 

for  all  p; 
(v)   if  h  is  real  and  (h,L  h)  ■  0  then  h  -  Oj 
' (vi;  f  is  purely  real  for  all  p. 
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With  assumptions  (iv)  and  (vi)  it  is  easy  to  see  that  all  g  and 
g'^  are  real.  Hence  if  (16)  holds,  assumption  (v)  shows  that  G    »  ^m+I  * 
In  particular,  this  applies  if  M  ■  -!• 

Consequently  we  have  shown: 

Theorem  3  The  sufficient  conditions  for  (f,G),  G  satisfying  the 
reciprocity  relation,  to  give  correctly  the  first  2M+2  terms  of 
the  power  series  expansion  of  (f,g)  are 


g"*  -  gj^  (m  <  M) 


when  assumptions  (i)  to  (iii)  are  satisfied.  If,  in  addition, 
assumptions  (iv)  to  (vi)  hold,  the  conditions  are  also  necessary. 

The  corresponding  theorem  for  the  more  general  case  is  not  so  powerful.  It  is 

•  I.        t 

Theorem  h    The  sufficient  conditions  for  (f  ,G)  or  (f,G  ),  G  and  G 

satisfying  the  reciprocity  relations,  to  give  correctly  the  first 

2M+2  terms  of  the  power  series  expansion  of  A  are 

g"'.  g^,  g'"-  g^  (iii<M) 

when  assumptions  (i)  to  (iii)  hold  with  appropriate  additicn  to  (i) 

I      t 
for  G  and  g  . 

The  proof  is  similar  to  that  for  Theorem  3  and  will  not  be  given. 

U,  The  low -frequency  behavior 

In  this  section  we  shall  consider  the  case  f  »  f  irtien  assumptions  (i) 
to  (vi)  hold*  These  assumptions  are  specially  applicable  to  the  low-frequency 
behavior  and  we  shall  make  particviiar  reference  to  the  process  employed  by 

r2i 

Levine  and  Schwinger '-  -* . 

n 
Assume  that  G 


■  3^  c  G   where  the  G  are  real,  independent  of  a 
m"0 


and  defined  in  a  similar  way  to  those  of  Theorem  1.  The  constants  c_  are 
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determined  by  equations  (9).  One  hopes  to  obtain  an  approximation  which  is 
reasonable  over  the  whole  frequency  band  by  a  suitable  choice  of  0  •  By  m%ans 
of  Theorem  3  we  can  say  what  this  choice  must  be  for  the  low- frequency  end. 
While  we  cannot  guarantee  from  this  work  that  the  choice  will  also  be  appro- 
priate for  the  high  frequency  end,  we  can  say  that,  at  any  rate,  it  is  not 
precluded* 

Suppose  that  we  desire  that  this  approximation  G  gives  the  first 
2n  +  2N  ♦  2  terms  of  (f,g)  correctly  at  low  frequencies,  N  being  an  integer. 
If,  in  particular,  n  ■  0  we  want  to  make  N  as  large  as  possible*  Now,  when 
the  first  2n  +  2N  +  2  terms  of  (f,g)  are  correct,  it  follows  from  Theorem  3 
that  G  »g^(m<n  +  N).  But  G  can  be  expressed  in  terms  of  G  ,G  ,.,,,G  • 
Hence  the  first  2n  +  2N  +  2  terms  of  (f,g)  can  be  correct  only  if  G  ,G-, o...G 

are  linear  combinations  of  go»^l»***'^n+N  °^»  alternatively,  if  go^Sp  •••jg^+jj 

are  linear  combinations  of  G  ,G,,...,G  . 

o'  1'   '  n 

If,  in  addition,  we  require  the  first  2n  +  2N  +  2  terms  to  be  corroct 
for  all  n  then,  on  taking  n  «  0,  g  ,g^,..»,g  must  be  constant  multiples  of  G  • 
Taking  n  ■  1,  we  deduce  that  g„^  must  be  a  linear  combination  of  G  and  G..  From 
n  "  2,  g«.^  must  be  a  linear  combination  of  G  ,  G_  and  G_  and  the  process  con- 
tinues. 

Thus  the  maximvm  number  of  terms  which  can  be  obtained  by  a  one- element 
approximation  to  G  (n  »  0)  is  determined  by  the  first  g  which  is  not  a  constant 
multiple  of  g  . 

Consider  now  the  application  of  these  results  to  the  problem  considered 
by  Levine  and  Schwinger.  In  that  problem  f-  ■  0,  L  a  0,  so  that  g,  ■  0,  but 
g-  is  not  a  constant  multiple  of  g  ,  Hence  N  ■  1  and  the  maximum  number  of  terms 
of   (f,g)  to  be  obtained  from  the  approximation  G  =  c  G  is  U  and  can  be  obtained 
only  by  choosing  G  ■  g  •  From  G  «  c  G  •♦•  c^G.  we  can  obtain  6  terms  provided 

OO  0011 
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that  G,  is  a  linear  combination  of  g  and  g_.  Similar iy,  a  suitable  choice 
1  o     d 

of  Gp  ^-'■1  Sive  8  terms  of  (f,g)  and  so  on»  The  choice  of  G  made  by  Levine 
and  Schwinger  is  in  fact  the  one  dictated  by  the  above  accuracy  considerations. 
Of  course  it  is  not  necessary  to  make  their  particular  choice  of  G^  to  ensure 
low- frequency  acc\iracy;  we  can  choose  any  linear  combinations  of  g  and  gp  so 
long  as  both  are  included.  Similar  remarks  apply  to  any  G  and,  in  general, 
we  woiild  choose  those  that  are  most  convenient  for  dealing  with  equations  (9). 

Finally,  there  is  one  thing  to  be  careful  of  in  the  use  of  Theorem  3» 
The  statement  of  accuracy  in  this  theorem  depends  upon  L  and  therefore  upon 
the  accuracy  with  which  L  is  approximated.  Thus,  when  the  inner  products  in 
(9)  have  to  be  evaluated  approximately  (as  often  happens),  care  must  be  taken 
that  this  evaluation  does  not  vitiate  Theorem  3»  Only  if  the  constants  in  (9) 

a    j  at  the  low -frequency  end  is  Theorem  3  available  . 

5.  Bounds  for  the  error 

Although  the  analysis  of  the  preceding  section  indicates  how  the  trial 
function  must  be  chosen,  it  does  not  give  a  precise  estimate  of  the  error  nor 
does  it  state  the  range  of  a  for  which  the  approximation  is  good.  In  practice, 
it  may  well  be  more  desirable  to  use  an  approximate  field  for  which  it  is  possible 
to  say  that  (f,g)  is  estimated  correctly  to  (say)  10°/o  for  a  certain  range  of 
wavelength  than  to  insist  that  the  field  satisfies  the  reciprocity  theorem  and 
gives  (say)  1  /o  accuracy  over  an  unknown  range.  We  shall  now  show  how  an  estim- 
ate of  the  error  can  be  made  in  certain  circumstances. 

Let  Y  ■  g  -G  where  G  is  some  approximation  to  g»  Then 

(f,Y)  -  (l€,Y)  •  (g,LY) 
or 
(17)  (f,g)  -  (f,G)  ♦  (G,LY)  +  (y,Ly)  . 
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Since 

(18)  LY  -  f  -  LG 

the  second  term  may  also  be  written  as  (G,  f-IG).  It  is  easy  to  see  from 
(17)  that  2(f,g)  -  (g,l€)  is  a  stationary  expression  for  (f,g)« 

If  G  satisfies  the  reciprocity  theorem  the  second  tenn  in  (1?)  is  zero 
and  the  error  made  by  the  approximation  is  (y>Ly)»  This  result  sometimes  en- 
ables us  to  state  the  sign  of  the  error.  For  example,  if  all  quantities  are 
real  and  L  is  positive  definite,  we  see  that  (f,G)  is  a  lower  bound  to  (f,g)« 

Even  if  G  does  not  satisfy  the  reciprocity  theorem  the  only  terra  which 
cannot  be  calculated  in  (17)  is  the  third  because  the  second  can  always  be  eval- 
uated by  means  of  (18).  An  estimate  of  the  goodness  of  the  approximation  there- 
fore requires  only  an  estimate  of  the  third  term  in  (17) J  this  estimate  depends 
essentially  on  being  able  to  deduce  properties  of  Y  from  a  knowledge  of  LY« 
There  does  not  appear  to  be  any  method  which  we  can  say  will  be  generally  applic- 
able because  the  operator  and  functions  occurring  are  usually  complex  but  we 
shall  give  two  methods  vrtiich  are  avsdlable  in  suitable  circumstances. 

Firstly,  suppose  that  it  is  possible  to  say 

1 2    /,   iu*\/i-  ^*\ 


i\»\)r  <     (h^,h^)(h2,h2) 


\rtiere  the  star  indicates  a  complex  conjugate  and  suppose,   fiirther,   that  there 
is  a  (i  >  0  such  that 

(LY,L*Y*)     >     ti(Y,Y*) 
for  all  Y  which  can  arise.     Then 


|(Y,LY)|     <     -i-     (Ly,lV)  . 


-  lU  - 


Sinc«  [L   and  LY  ar  both  known  this  inequality  gives  an  estimate  of  the  error. 
We  may  also  determine  the  accuracy  of  the  other  approximations.  For 
(f',  g)  -  (f',  G)  ♦  (f',  y)  and 


(19)  |(f  ,Y)r  1  ^  (f  ,f  )(LY,L  Y  )  } 

all  quantities  are  known  on  the  right  hand  side  of  (19 )♦ 

A  second  method,  related  to  that  used  by  Kato*-  -•  in  problems  of 
atomic  collisions,  depends  upon  adjusting  the  incident  field  until  a  pre- 
scribed scattered  amplitude  is  obtained.  The  scattered  amplitude  correspond- 
ing to  the  original  incident  field  can  then  be  deduced  immediately.  Let  v  be 
a  soFp'lex  constan"  such  that,  if 

(20) 
(21) 


V  f  - 

Lg°, 

(f. 

,g°)  -  1  . 

'  «  vg 

we 

see 

that 

(f,g)  ' 

1 
'  V  • 

Consider 

now 

Q  chosen  so 

that 

(22)  (f,G)  -  1  , 

Let  G°  "  V  f  -  (L-1)  G  -  V  f  -  Mj  where  v  is  a  complex  constant,  G  is  con- 
sidered as  an  approximation  to  g  •  Then 


(g°,  G  -  G°)  -  (g°,  LG  -  V  f )  -  (l€°,  G)  -  V 

t 

■   V  -  V 

on  account  of  (20,  (21)  and  (22),  Hence 

(23)  V  -  v'  +  (G°,  G  -  G°)  ♦  (g°  -  G°,  G  -  G°)  , 
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We  can  regard  the  first  two  terms  as  giving  an  approximation  to  v,  the  third 
being  the  error. 

To  estimate  this  error,  we  introduce  the  equation 

(2U)  (l  +  \)Mj^  +  0  -  j\.f 

where  (1  +  X)(f,0)  ■  0  and  _A. is  a  non-zero  constant  to  be  chosen  later.  We 
assume  that  non-zero  0  can  be  obtained  only  for  a  discrete  set  of  real  \  of 
which  X  is  a  typical  one,  0  being  the  corresponding  solution*  Since 

1  +  X 

(0  -^  )  -  -(X  +  i)(0  ,  M  g(  )  -  -(X  +  i)(M0  0  )  -  r-r-^    (0  »   '^  - 


n 


we  choose  _/V  for  each  X  so  that 

n 


(0,0)  -  5 
^'^m'  '^n      nra 


where  6   ■  0  (n  /  ra),  «  l(n  -  m).  It  will  be  assumed  that  the  0  form  a 

complete  set. 

Define  the  constants  a     and  b     by 

n     n  "^ 

(G  -  G°,  0  )  -  a 
^     '  '^n    n 


and 


Then 


(g°  -  0°.  0„)  .  b„  . 


b  -  (vf  -  M  g°  -  v'f  +  M  G,  0  ) 


(M  (G  -  g°),  0^) 


because  of  (2U)j  thus 
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b^    -     (G  -  g°,   M  0^) 


on  account  of   (21)  and   (22).     Also 


'TTT    ^''  ^°>  ^n) 


a^     -      (IC  -  l€°,   0^) 


(G  -  g°,    L  i?^) 


(G  -  g°,   0 J 


^*\  '      ^ 


-    -X„  b     . 
n    n 


Therefore 


(g°  -  G^  G  -  G°)    -   L  ^  ^  -  -E  r-  4 


and 


2 


(G-G,G.G)     -     Yi     < 


• 


If,   now,  we  can  say  that  a     is  real  for  all  n,   it  follows  that 

-  ^  (G  -  G°,   G  -  G°)  <   (g°  -  G°,   G  -  G°)  <  -  i  (G  -  G%  G  -  G°) 

where  X^  is  the  smallest  positive  X  and  X__  is  the  negative  X  of  smallest  ab- 
solute value.  (Note  that  X  cannot  equal  zero  for  any  n. )  Since  G  is  known 

t 

in  terms  of  v  and  G  we  have  obtained  an  estimate  for  the  error  in  using  the 

first  two  terms  of  (23)  as  an  approximation  to  v.  The  accuracy  of  G  as  an 

o 

approximation  to  g  may  also  be  obtained  from 
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O    „0    O    „0^  T—   ,2 


<  max 


This  enables  us  to  determine  the  error  in  quantities  such  as  (f  ,  g  ). 

In  practice  this  method  may  not  be  at  ail  easy  to  apply  because  of 

the  many  restrictions.  We  have  first  to  verify,  in  any  particular  problem, 

that  the  fJ  s  satisfy  the  conditions  imposed  and  then  that  G  and  v  can  be 

chosen  so  that  a  is  real,  (The  latter  condition  is,  at  any  rate,  complied 

with  when  f  and  L  are  real.)  Then  X^  and  X__  (or  bounds  for  them)  must  be 

detennined.  Fortunately,  great  accuracy  is  not  required  here  because  we  make 

(G  -  G  ,  G  -  G  )  as  small  as  possible  by  taking  the  best  choice  of  v  among 

those  values  which  keep  a  real, 

n 
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